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lO ' Abstract 
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^ I The rare decay ^ 77 is a promising detection channel for an intermediate mass 

p ■ Higgs boson. We compute its two- loop Ola'^GF'mf) correction in the standard model 

CL|' and find that the relative correction to the decay rate runs between 0.7% and 0.5% for 

^ • Mh = 80 — 150 GeV. The analogous correction to the amplitude for gg H is recovered 
as a special case. The generalization of our result to other models is also briefly indicated. 
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The standard model ( SM ) of the electroweak interactions has proved very suc- 
cessful in the description of electroweak phenomenology. Yet, there is one particle, the 
Higgs boson, predicted by the SM that has evaded detection up to now. The Higgs boson 
is an essential ingredient of the SM. It provides mass for the W, Z bosons and fermions 
through the Higgs mechanism. Its discovery will thus be crucial for confirmation of the 
Higgs sector in the SM. 

The search for the Higgs boson is difficult however. Theoretically this is basically 
because the relevant parameters such as masses, couplings are essentially free in the SM. 
Although there are theoretical considerations that can more or less constrain them our 
knowledge about these parameters is mainly from the failure search of experiments. For 
example, experiments at LEP I and SLC have ruled out the Higgs mass range < 63.5 
GeV at the 95% confidence level. It is expected that LEP II will extend the range up to 
80 GeV. Beyond this we have to appeal to the next generation of hadron colliders. The 
production mechanisms for the Higgs boson at hadron colliders have been extensively 
studied in the literature 0. Generally the gluon fusion mechanism [Q] dominates for 
a Higgs mass up to 700 GeV. Above this range the dominant mechanism is through 
W, Z scattering subprocesses. The detection mechanisms for Higgs boson in the range 
Mh ~ 140 - 800 GeV have also been widely studied ( see, for example, Refs.[|[@|[|]). 
Less extensively studied is the so-called intermediate mass Higgs boson in the region 
Mh ~ 80 — 140 GeV 0. It is difficult to detect a Higgs boson in this region because the 
dominant decay mode H bb is badly buried in the enormous QCD jet background]^ so 
that one has to use rare decays. A favorite decay mode is H ^ 77, though an excellent 
energy resolution is still required to discriminate signals from the background produced 
from qq{gg) 77 and fake 7's from 7r° decays^]. Since the signal of if 77 is so small 
and may be overwhelmed by the background it is important to predict its decay rate as 
precisely as possible. 
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At lowest order the amplitude for the decay — 77 receives contributions from 
charged fermion and W boson loops^]. Among fermions in the SM, only the top quark 
is of practical importance in fermion loops0. The QCD correction to the top loop was 
considered in Ref. and found to be well under control. The 0{a'^GFiTi'f) correction 



from internal exchange of the Higgs boson in the top loop was computed in Ref. |T^. In 
this paper, we will complete this by including contributions from unphysical Goldstone 
bosons as well. The gluon fusion process gg ^ H is quite similar to the decay H 77. 



The QCD correction to the fusion was shown to be very large |Tl| and this would make 
the 0{Gpm'^) corrections relatively more important for iJ ^ 77 than for gg —>■ H. 
The Ola'^Gprnf) correction to the fusion has been computed by the use of the low energy 
theorem for the trace of the energy-momentum tensor and found to be small[0. Although 
the 0{alGFni1) correction for gg ^ H may be reproduced from that for i7 — 77 by 
assigning equal charges to the top and bottom quarks, the latter cannot be obtained 
simply from the former as shown below. 

The lowest order contribution to the amplitude for if ^ 77 is 

O = ^el{l)e:{2){g^^k, ■ k, - k^k'(), 
Zttv ^ 

where ki 2 and e(l, 2) are the momenta and polarization vectors of the two photons, Qt^b 
are the electric charges of the top and bottom quarks, is the number of colors, and 
V = 2-^/^Gp^''^ = 246 GeV. Ft and Fw are fuctions of r]t = Amj/Ml and i]^ = 4M^/Mfj 
respectively, given in Refs. To obtain the desired O^a^Gp'ml) correction, we 

work in the heavy top limit. This is a good approximation for an intermediate mass Higgs 
boson. In explicitly renormalizable gauges, the leading term in this limit is provided 
by internal exchange of the Higgs and unphysical Goldstone bosons with minimal number 
of contact interactions amongst themselves. This is because the introduction of a contact 
scalar vertex will bring in a factor of Mfj = 2ki ■ k2 that cannot be cancelled by the 
possible logarithmic infrared behaviour of diagrams in the limit mf ^ Mfj, and thus 



will not contribute to the leading term. Furthermore, the non-leading terms are ^w,z- 
dependent, and this dependence is cancelled only when contributions from W, Z bosons 
are included. With these considerations, the relevant interaction Lagrangian is 

V V V iZ) 

+te{Qt - Qb)A^{<j)-dy+ - + e\Qt - g,)2A^A^0+0-, 

where H, cfp'^ are the Higgs and unphysical Goldstone bosons respectively, and ^4^ the 

photon field. We have ignored the small bottom mass and quark mixing. 

The two-loop diagrams that may contain the desired correction are obtained by at- 
taching in all possible ways the two photon lines onto the diagram shown in Fig. 1. They 
are shown in Fig. 2. These diagrams are classified into three groups. The first group 
corresponds to the insertion of the following one-loop elements in the one-loop top or 
bottom diagram for H —>■ ■yj: the top self-energy (denoted asW ), the Hti or Hbb vertex 
(Vi), and the A^tt vertex ( V2 and V3 ) . The second group consists of the insertion of the 
one-loop H(l)'^(/)~ vertex in the one- loop 0^ diagrams for — > 77 ( V4 ). The third group 
is the remaining overlapping diagram { R). Note that diagrams Vs,4, and R are present 
only when the exchanged virtual scalars are the charged one, 0=*=. It is these diagrams that 
will produce the difference between the 0{GFmf) corrections to if ^ 77 and gg H. 

The contributions from the i?— , and 0^— exchanged diagrams are respectively 
gauge invariant and can be calculated separately. In addition, parity is violated by the 
(f)^th vertex so the 0^— exchanged diagrams will generally induce a new parity- violating 
Lorentz structure, O^^ = ien^p„k{k2. Since the one-loop amplitude is symmetric and O^,^ 
antisymmetric with respect to // and the parity-violating term will not contribute to 
the total decay rate to the order considered here and may be safely ignored. In other 
words, there are no ambiguities associated with the 75 problem though we are computing 
two-loop diagrams. 

Now we describe the renormalization procedures used to define the physical parameters 
and renormalization constants. We use the dimensional regularization to regulate ultravi- 
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olate divergences and work in the on-shell renormalization scheme. For our purpose, the 
mass of the unphysical Goldstone bosons can be set to zero from the very beginning, but 
with one caveat as explained below. 

1. Mass and wavefunction renormalization of the top { W ) - The counterterms for 
the top mass ( Srrit ) and wavefunction renormalization constant ( SZt ) are determined 
by requiring that rrit be the pole of the one-loop corrected propagator of the top and 
that the residue of the propagator at its pole be unity. Note that in the case of 0^ 
exchange there are wavefunction renormalization constants for the left- and right-handed 
parts respectively, ^ 1 + 5Z^, Z^ ^ 1 + 5Z^. 

2. Renormalization of the Hti and Hbb vertices ( Vi ) - In the SM, this is not indepen- 
dent but is related to the renormalization of the t, b self-energy. The Feynman rule for the 
counterterm of the Htt vertex induced by H— or exchange is —imt/v{5mt/mt-\-5Zt). 
For (jy^ exchange there is no counterterm for the Hbb vertex in the limit of zero bottom 
mass, but there is a counterterm for the Hti vertex, —imt/v{8mt/mt -\- 1/2{SZ^ + Zf)) 
though its bare one-loop diagram does not exist in this hmit. Finally there is a global 
counterterm for the Htt vertex, [—imt/v){5v/v) where 5v is the VEV counterterm of the 
Higgs field induced by the heavy top. 

3. Renormahzation of the A^jti vertices ( 1^2,3 ) - The electric charge is defined as 
usual in the Thomson limit. We have checked that the f/(l)c.m. Ward identity is satisfied 
especially in the case of 0^— exchange where this is not self-evident. 

4. Renormalization of the H(f)'^(f)~ vertex ( V4 ) - The computation of diagrams V4 
requires special care. Let us first discuss the renormalization of the vertex. The countert- 
erm for the vertex is —i{M]j/v){SZx — Sv/v)ii'^, where SZx is the renormalization constant 
for the scalar self-couphng A, 5Zx = 5Zh + {l/M%) {5M% - 5Mj±) + 25v/v. The mass 
and wavefunction renormalization of the Higgs boson ( SMjj, 6Zh ) is done as for the top 
so that there is no renormalization factor for the external Higgs boson. The counterterm 
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for the mass of 0^ ( 5M|± ) is determined by the condition of tadpole cancellation at the 
one-loop level. The masslessness of cfP'^ is then automatically preserved by the Goldstone 
theorem at the same level. Since the counterterm for V4 is to be inserted into a finite 
one-loop amplitude and only terms up to 0{rn^) are required, we obtain 



5v 7 2 1 

where — = -NA- — ) has been mserted. 
V 6 A.t:v 

It is straightforward to calculate the contributions from H— or 0°-exchanged diagrams, 
i.e., diagrams W and Vi 2. The bare H-exchanged diagrams sum to a finite, gauge invariant 
form. In the heavy top limit, it is given by 

^At:^.)=^ON^{^f[-m (4) 
The counterterm diagrams as a whole are also finite and gauge invariant, 

where 5m\^^ is the counterterm for the top mass induced by H-exchange, and i^^"'""^ 
must be computed in n-dimensions, i.e., 

Ft--\n^+eW^r. (6) 

Note that the wavefunction renormalization constant for the top is cancelled. Similarly, 
we obtain. 

It is much more difficult to calculate the contributions from 0^-exchanged diagrams. 
Besides Vi^2,3 and i?, they involve V4 that are seemingly infrared divergent in the limit 
of massless 0^. Although they are actually not infrared divergent, individual diagrams 
do contain terms that are non-analytic in = 2ki ■ k2 ( Higgs boson momentum squared 



), like fif^i^ln"-^ {n = 1,2), gau^\n^, ^^2ith!L^ hfj,hu P i^^^i^P^ jg 

mf mf mf mf mf 

a non-trivial check of our calculation that these terms are cancelled in the sum of two 
diagrams in V4 so that they will not spoil our low energy expansion in the heavy top limit. 
The imaginary part is also cancelled in the sum. This is reminiscent of the observation 
that the one-loop 0^-exchanged amplitude for iJ — 77 does not contain an imaginary 
part in the limit of massless 0"*^. The sum of all bare 0^-exchanged diagrams and the 
counterterms for V4 is 

^•Kn^Ze) = ^ON^{^n4Q,Q, - ^Q?) + 5{Q, - Q,)'], (8) 

where the second term is contributed by diagrams V4 and countermterms, and the first 
term is contributed by other bare diagrams. We note that diagrams V4 actually contain 

an additional term proportional to (Qt — Qb)^guu\—^ — 4:A2e ^1, which is exactly 

cancelled by other diagrams. The remaining counterterms also sum to a gauge invariant 
form, 

= Sraf^^Ul-"-' = ^ON.{^fl|Q^ (9) 

Notice that the wavefunction renormalization constants and 6Z^ are again cancelled 
in the sum. This is because the absence of a counterterm for the Hbb vertex is compensated 
for by the presence of a counterterm for the Hti vertex. The same counterterm also makes 
the drrif^ •* part just as simple as in the case of H or (f)^ exchange. 

Finally, there is a contribution from the counterterm 6v/v for the Htt vertex, 

ML-^°°p = = ^ON^{^n-^-^N,Q^ (10) 

Three different methods are employed to compute two loop diagrams. (1) After loop 
integration we expand Feynman parameter integrals in the heavy top limit to obtain the 
Qnu and k2nkiy terms. (2) We slightly modify the Hoogeveen's method on expansion in 
the large mass limit |]l3i. We directly expand top propagators to the desired order without 



shifting 7 matrices from denominators to numerators beforehand, and then use algebraic 
identities to further reduce their products. This simphfies algebra considerably. We use 
this to get the k2^kii, terms. (3) After loop integration we use numerical analysis to 
approach the heavy top limit. Both (7^1, and k2^kiij terms are computed. In the case 
of H-exchange all three methods lead to an identical result. For 0°-exchanged diagrams 
we use the first two methods and indeed obtain the same numbers. It is complicated to 
apply the method (2) to 0''=-exchanged diagrams due to the infrared behaviour associated 
with masslessness of the bottom and 0^, so only the method (1) is used. But even so, 
we still have nontrivial checks as mentioned above: cancellation of non-analytic terms 
in diagrams V4, and cancellation of divergent, non-gauge-invariant terms proportional to 
QnuiQt — QbY between V4 and other diagrams. To summarize, the complete O^a^Gp'Tnl) 
contribution to the decay amplitude is 

The amplitude for the fusion gg H is recovered by setting Qt = Qb and changing 
coupling factors appropriately, 
i-Aiga 9b H) 

= ^e,(l)e.(2)(,-fc, . k, - A;,X)tr(^^)(- J)[l + {^ns + ^iVj] (^3) 

= ^e,{l)e^mg'^^k, ■ k, - k!^k-(M^^){--)[l + ^^3^]^ 
27rv ^ 2 1; V 2 2 3^^ 327r2 ^' 

which coincides with the result of Ref. [O. i 



Including 0{a^GF^1) and QCD corrections, the decay rate for if 77 is 



r = r,_,„„|i + ^%^iB + H;!ic]. (13) 



where Fi.ioop is the lowest one loop contribution, C was computed in Ref. and 

14 

B = N,[AQtQ, + 5(Qt - Q,f - —N,Ql]l[N,QlFt + {Q^ - Q.fFw]. (14) 

y 



^We thank A. Djouadi for pointing out an error in our comments on Ref. |12[ in the original version 
of this paper. 
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For numerical analysis we use rrit = 176 GeV, then the 0{a^GFm^) correction to the 
decay rate runs between 0.7% and 0.5% for Mh = 80 — 150 GeV, which is roughly one 
half of the corresponding QCD correction^ if ~ 0.1 is used. 

The results reported here may be employed to incorporate contributions from exchange 
of extra scalars in models with an extended Higgs sector. For example, in the two Higgs 
doublet model [1^ the contributions from all of the four physical scalars may be incor- 



porated by multiplying appropriate factors coming from vertices. But it is suspect that 
the heavy top limit remains to be a good approximation as these particles are generally 
much heavier than the top quark. 

We thank K.-T. Chao and Y.-P. Kuang for discussions. 
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Figure Captions 

Fig. 1 The two-loop diagrams are obtained by attaching in all possible ways the two 
photon lines onto the diagram shown here. Solid and dashed lines represent fermion and 
scalar fields respectively. 

Fig. 2 The two-loop diagrams to be computed here. Wavy lines represent photon fields. 
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